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SUMMAR 


Heat  transfer  and  pressure  distributions  have  been 
made  on  five  vehicle  nose  shapes  in  the  VKI  Longshot  free  piston 
wind  tunnel  at  H«15  and  20.  The  flov  parameters  achieved  closely 
simulate  aerodynamic  re-entry  conditions.  The  models  used  in 
the  study  were  a  50°-8°  biconic  configuration,  vith  sharp-  and 
blunt-nosed  versions  and  a  sharp-nosed  version  vith  a  machined 
roughness,  and  a  hemisphere  configuration  vith  a  smooth  and  a 
sand-blast  rough  surface.  These  configurations  resemble  the 
stable  shapes  of  turbulent  and  laminar  ablating  nose  cones. 
Preliminary  comparisons  have  been  made  vith  the  stai&e-of-the  art 
engineering  predictions.  Some  discussion  about  the  state  of 
the  boundary  layer  is  included. 
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1.  INTRODUCTION 

1 

Many  investigations  of  the  heat  transfer  and  pressure 
distributions  on  re-entry  nose  shapes  have  been  made  in  the  p&st. 
However  little  experimental  data  exists  on  models  tested  in  truly 
re-entry  simulated  conditions  of  high  Mach  number  and  high  Reynolds 
number  conditions.  Measurements  of  heat  transfer  on  smooth  ahd 
rough  hemispherical  and  conical  bodies  will  aid  predictions  of 
the  performance  of  tne  ablation  method  of  heat  protection.* 

A  laminar  ablating  nose  takes  on  a  hemispherical  stable  shape. 

In  the  turbulent  case  a  biconic  shape  with  roughly  a  50®  fore¬ 
body  with  an  8°  degree  skirt  is  formed  after  a  long  period  (Ref.  l). 

Some  of  the  outstanding  problems  necessary  to  predict 
ablati'on  recession  rates,  for  instance,  are  involved  with  not 
only  verifying  the  theories  used  for  heat  transfer  prediction  but 
also  knowing  the  location  of  the  transition  of  the  lamina^  boun¬ 
dary  layer  to  a  turbulent  boundary  layer  on  bqth  smooth  and  rough 

i 

blunt  forebodier.  1  > 

The  current  state  of  the  art  techniques  for  predicting 
laminar  heat  transfer  rates  over  nose  shapes  (R6f.  1)  involve 

i 

the  use  of  the  laminar  stagnation  puint  correlation  of  Fay  and 
Riddell  (Ref.  2)  with  a  local  similarity  solution  (Ref.  3).  For 
turbulent  flows  then  the  reference  enthalpy  method  of  Eckert 
(Ref.  k)  or  Spalding-Chi  (Ref.  5)  for  flat  plates  is  used  or  a  . 
calculation  based  on  the  Vaglio-Laurin  (Ref.  6)  correlation 
which  accounts  for  pressure  gradients. 

i 

The  models  used  in  experimental  study  are  cf  two  basic 

i 

types,  a  hemisphere  and  a  50°-8°  biconic  configuration.  Both  a 
smooth  and  an  intermediate  (sand-blasted)  roughness  surface  was 
used  in  the  former  case.  For  the  biconic  shape  both  sharp-  and 
blunt-nosed  smooth  models  were  tested  as  well  as  a  sharp-nosed 
version  with  a  very  deep  waffle-iron-like  roughness  machined 

.  I 

into  its  forebody.  Heat  transfer  and  pressure  distributions  were 

I- 
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measured  on  these  models  in  the  M«20  f lov  of  the  Y.K.I,  LongBhot 

•  , 

* . unel,  The  aerodynamic  conditions  in  this  vind  tuntael  nearly 
simulate  the  conditions  actually  obtained  on  a  re-entering  missile 
during  its  most  critical  maximum  heating  and  deceleration  phase* 

•  i  '  • 

The  temperatures  at  the  outside  edges  of  the  bopndary 

i  .  1 

layers  formed  on  these  blunt'  models  are  very  high,  of  the  order 
of' 2000-2500®K,  compared  to  the  vail  temperature  at  ambient  con¬ 
ditions*  These  aerodynamic  conditions  iare  beyond  those  previously 
compared  with  the  above-mentioned  predictions  for  heat  transfer. 

The  gas  in  the  boundary  layer  ds  in  a  vdbrationally  excited  state* 
Because  of  the  h'igh  temperature,  the  Reynolds  numbers  on  the  sur¬ 
face  are  lov,  despite  the  high  Reynolds  number  of  the  stream. 

Under  such  advanced  and  realistic  conditions  achieved  in  these 
tests,  it  vill  be  possible  to  check  the  ability  of  the  current 
theories  to  predict  the  heat  transfer  rate  on  nose  cones  ablating 
1  in  a  high  enthalpy  stream.  , 

!  i 

i 

2,  EXPBRlkESTAL  APPARATUS. 

2,1  Models 

i  ' 

Five  instruaented  models  fabricated  and  supplied  by 
AYCO  weW  tested  in  the  Lcngshot  floe1.  Theae  models  vere  deri¬ 
vatives  of  tvo  basic  shapes,  one  a  biconic  and  the  other  a  hemi- 

t  ,  i 

spherical  model.  The  fc iconic  model  had  a  50°  half-angled  conical 
forebody  of  4  in  surface  length  and  an  8C  half-angled  afterbody 
leading  ,to  a  7  in  diameter  base.  The  hemisphere  model  had  a 
diameter  of  7  in?  ! 

l 

The  biconic  shape  vas  tested  in  three  form*,  k  sharp- 
nosod  model  vitb  a  saeoth  vail  (called  model  A  in  this  report), 
a  sharp-nosed  model  vith  the  forebody  roughened  vith  a  vaffle- 

*  i 

iron  pattern  machined  deeply  int®  the  surface  (B)  and  a  blunt- 
nosed  model  (C)  vith  a  nose  radius  of  0,75  in  and  a  smooth  sur- 

i  ,  * 

face.  The  hemisphere  vas  tested  in  tvo  fbr^s.  One  a  smooth  model  (D) 


“  J  - 


a 


and  the  other  a  model  roughened  with  a  sand-blasted  finish  (E). 
Photographs  of  models  B,  C,  D  and  £  are  shown  in  Pig#  1. 
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The  models  were  fitted  with  pressure  taps  and  heat 
transfer  gauges,  whose  positions  are  shown  in  Pig.  2  and  listed 
in  Table  1.  For  each  model  the  meridian  containing  the  pressure 
taps  was  at  180  degrees  from  the  meridian  containing  the  heat 
flux  gauges.  Because  instrumentation  limitations  only  10  heat 
transfer  measurements  and  7  pressure  measurements  could  be  made 
in  one  Longshot  test. 


2.2  Beat  transfer  instrumentation 

The  flux  gauges  were  fabricated  of  O.OC^  j.n  thick  copper 
disc  bonded  to  an  insulating  holder  as  shown  in  Fig.  3.  A  chromel- 
alumel  thermocouple  junction  approximately  0.0013  inches  in  dia¬ 
meter  vas  formed  by  spot  welding  the  0.001  inch  diameter  thermo¬ 
couples  to  the  backface  of  the  copper  disc.  All  gauges  were  con¬ 
toured  to  the  local  body  configuration.  The  discs  had  a  diameter 
of  0.125  inches  except  for  the  rough  biconic  model  ( B )  which  had 
gauges  with  disc  diameters  of  0.110  inches. 

Calibration  of  typical  heat  flux  assemblies  were  ob¬ 
tained  in  the  AEX>e  radiant  heat  flux  calibration  facility.  The 
purpose  of  the  calibration  was  to  obtain,  for  a  given  set  of 
standard  material  properties,  the  gauge  assembly  effective  thick¬ 
ness  to  be  used  in  the  heat  flux  data  reduction  analysis.  Of  six 
gauges  tested  the  average  nominal  percent  deviation  is  a  reduc¬ 
tion  of  5$  in  thickness,  although  the  maximum  deviation  varied 
fros  to  -12.6 $  in  thickness. 

? >r  determining  the  calibration  factor  of  a  gauge,  the 
following  material  properties,  supplied  by  AVCO,  of  ©wygen-free 
high  conductivity  copper  and  cbromel-aluael  thermocouples  were 
used.  Ths  values  were  averaged  over  a  temperature  range  from 
295° X  to  3£5°K, 


-  u  - 

Density  “  558  lb/ft3 

Specific  heat  •  0,915  B,Th.U/lb/°F 
dT/dV  -  kk,2°F/mV 


Using  Equation  B,2  (Appendix  B)  the  calibration  cons¬ 
tant  for  a  calorimeter  vith  a  disc  of  0,00b  inch  thick  vill  be 

0,752  /  'SiL-.  .taking  the  5S  reduction  in  thickness  into 

_ .  o  sec 

ft^sec 

account  the  calibration  constant  for  aO.OQbinch  thick  gauge  used 
in  the  data  reduction  programme  (see  Appendix  B)  vas  0,715,  Gauges 
with  discs  of  differing  measured  thickness  were  corrected  propor¬ 
tionally  to  that  dimension. 
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Output  signals  from  the  calorimeter  gauges  were  generally 
recorded  on  Tektronix  oscilloscopes  using  Polaroid  cameras.  Pour 
channels  used  some  pre-amplification  using  D,C,  amplifiers  built 
at  V.K.I,  Tvo  chancels  used  in  tests  on  models  A  and  B  were 
recorded  on  a  CEC  oscillograph  after  amplification  by  CEC  D.C. 
amplifiers.  All  this  instrumentation  vas  calibrated  at  frequent 
intervals. 

Early  tests  were  carried  out  in  Longshot  on  gauges  with 
stainless  steel  calorimeters  of  0,001  thickness  supplied  by  AVC0. 
These  tests  shoved  the  gauges  to  be  fragile.  The  heat  conduction 
losses  to  the  thermocouple  vires  vas  expected  to  cause  large 
corrections  to  the  theoretical  calibration.  These  gauges  were 
rejected  in  favour  of  copper  calorimeters  used  throughout  these  tests. 


2,3  Pressure  instrumentation 


v 


The  pressure  taps  mounted  on  the  model  had  internal 
diameters  of  in  and  channel  lengths  of  in.  The  length 

of  rubber  tubing  of  internal  diameter  in  connecting  the 

pressure  taps  to  the  transducer  vas  in.  The  internal  volume 

of  the  total  pressure  systems  vas  small  enough  to  give  a  pneu¬ 
matic  response  time  of  less  than  1  millisecond. 


2,4  Measurements  of  the  reservoir  conditions 


The  reservoir  pressure  vas  obtained  using  a  Kistler 
type  6201  gauge  mounted  in  the  vail.  The  signals  vere  processed 
and  recorded  using  a  Kistler  charge  amplifier  and  a  Tektronix 
oscilloscope  fitted  vith  a  Polaroid  camera. 

The  reservoir  temperature  vas  measured  vith  a  tungsten- 
25 %  rhenium  vs,  tungsten-3?  rhenium  thermocouple  aounted  in  the 
reservoir  vail.  The  vire  diameter  used  vas  0,010  in  vith  a  total 
length  of  approximately  0,5  in  giving  a  response  time  from  5  to 
10  milliseconds.  Conduction  loss  s  had  been  assessed  in  previous 
tests  by  using  thermocouples  vith  differing  t/d  ratios.  The  method 
of  data  reduction  is  described  in  Appendix  A. 

Typical  traces  from  all  the  above-mentioned  instrumen¬ 
tation  are  shovn  in  Fig.  5* 

2.5  Schlieren  system 

An  18  in  conventional,  single  pass  Toepler  schlieren 
system  equipped  vith  high  quality  optical  components  is  used. 

With  the  exception  one  24  in  plane  mirror  to  bend  the  l^ght 
90°  (due  to  the  vicinity  of  a  vail  near  the  test  section)  the 
light  beam  takes  a  Z  shaped  path.  A  single  spark  light  source 
vith  a  spark  duration  of  1  microsecond  is  used  to  record  the 
visualization  of  the  flov  on  either  a  Polaroid  or  a  sheet  film. 

The  magnification  of  the  system  is  approximately  0,25*  A  grid 
made  of  nylon  line  at  1  cm  intervals  is  placed  at  one  of  the 
observation  vindovs.  A  grid  is  thus  placed  on  each  schlieren 
photograph.  The  lines  are  slightly  out  of  focus,  since  the 
schlieren  system  is  focused  on  the  model. 
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2.6  The  Longshot  wind  tunnel 

The  teats  vere  carried  out  in  the  VKI  Longshot  free- 
piston  vind  tunnel  (Ref.  7).  A  photograph  of  this  facility  it 
shovn  in  Fig.  3.  Longshot  was  designed  by  Perry  and  his  asso¬ 
ciates  at  Republic  Aviation  specifically  for  the  attainment  of 
very  high  Reynolds  number  at  high  Mach  numbers.  It  differs  from 
a  conventional  gun  tunnel  in  that  a  heavy  piston  is  used  to 
compress  the  nitrogen  test  gas  to  very  high  pressures  end  tem¬ 
peratures.  The  test  gas  is  then  trapped  in  a  reservoir  at  peak 
conditions  by  the  closing  of  a  system  of  check  valves  as  the 
piston  rebounds.  The  flov  conditions  decay  monotonically  during 
10  to  20  milliseconds  running  time  as  the  nitrogen  trapped  in 
the  reservoir  flows  through  the  nosale  into  the  pre-evacuated 
open  jet  test  chamber.  The  maximum  supply  conditions  currently 
available  are  60,000  psi  at  2400°K.  Hossle  throat  inserts  giving 
nominal  Mach  numbers  of  15  and  20  at  the  nozxle  exit  of  l4  in 
diameter  vere  employed  in  these  tests. 

The  rough-hemisphere  model  is  shovn  mounted  in  its 
test  position  in  the  Longshot  working  section.  The  shroud  mounted 
around  the  model  is  placed  to  prevent  disturbances  arising  from 
reflections  of  the  flov  from  the  back  vail  of  the  dump  tank. 

The  shroud  has  been  proved  to  be  effective  in  preventing  distur¬ 
bances  in  the  test  section  occurring  before  40  milliseconds 
after  the  peak  conditions.  Such  disturbances  vere  mentioned  in 
Ref.  8.  The  channel  sections  attached  to  the  shroud  lead  to  the 
schlieren  observation  vindovs.  The  model  is  shovn  mounted  on  a 
sting  attached  to  a  quadrant  traversing  system. 
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3.  TEST  COHDITIOIS 

The  test  conditions  vere  computed  for  each  run  using 
the  Longshot  data  reduction  programme  (outlined  in  Appendix  A). 
Table  2  summarises  the  tests  carried  out  in  Longshot.  Measured 
reservoir  conditions  and  Pitot  pressures  at  peak  conditions  are 
given,  Also  tabulated  are  the  calculated  conditions  of  Mach 
number,  Reynolds  number,  dynamic  pressure  and  stagnation  point 
heat  transfer  on  a  7  in  diameter  hemisphere.  Three  basic  test 
section  conditions  vere  used  throughout.  These  vere  a  high  and 
a  lov  Reynolds  number  case  at  M*20,  and  a  high  Reynolds  number 
case  at  M*15.  The  lov  Reynolds  number  case,  however,  vas  repea¬ 
ted  at  several  slightly  differing  conditions.  This  vas  done  since 
this  conditions  vas  tried  for  the  first  time  during  the  tests, 
and  it  vas  vished  to  make  sure  that  the  lovest  Reynolds  number 
vas  achieved  vithout  condensation  occurring. 

It  is  seen  that  there  is  reasonable  repeatability  in 
the  data.  Any  discrepancies  are  as  much  due  to  the  compression 
cycle  (since  the  pressure  developed  is  very  sensitive  to  the 
barrel  pressure)  as  the  scatter  in  measurement.  The  temperature 
is  probably  the  least  accurately  measured.  Some  difficulties 
occurred  vith  the  Pitot  probe,  since  in  the  M*20  runs,  the  probe 
vas  located  very  close  to  the  norsle  turbulent  boundary  layer. 

In  some  cases  this  measurement  is  suspect. 

Fuller  details  of  the  test  section  conditions  are  given 
for  selected  runs  for  each  set  of  conditions  (Table  3).  The  sym¬ 
bols  and  units  are  explained  in  Appendix  A. 
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4.  RESULTS 


All  measurement*  of  the  peak  values  of  heat  transfer 
and  pressure  are  tabulated  in  Tables  4  to  7  and  Tables  8  to  1 1 
respectively.  Schlieren  photographs  of  all  tests  are  show  *  in 
Fig.  8.  The  results  from  each  model  configuration  vill  b&  dis¬ 
cussed  in  the  folloving  sub -sections . 

4.1  Pointed  smooth  biconic  model  (A) 


Examination  of  the  schlieren  photographs  of  Runs  204 

to  210  show  that  the  bov  shock  remains  attached  to  the  nose  for 

all  cases.  For  the  aero  angle-of-attack  cases  ,  the  shock  vave 

angle  is  approximately  7-5°  to  the  surface  which  agrees  to  within 
1  ° 

2  to  the  theoretical  results.  This  indicates  that  the  flow  is 


supersonic  behind  the  shock. 


At  the  0=410  cases  the  shock  shape  again  agrees  with 

1° 

the  theoretical  cone  tables  of  Jones  (Ref.  9)  tc  within  —  .  This 
is  in  spite  of  the  fact  that  the  tables  have  been  extrapolated 
to  such  an  extent  that  they  give  subsonic  conditions  on  the  wind¬ 
ward  surface.  On  the  leeward  surface,  She  flow  is  predicted  to 
remain  in  a  supersonic  condition.  (It  is  interesting  to  note  that 
for  0c>2O°e  the  shock  angle  to  the  surface  increases  with  increas 
ing  at  those  high  Mach  numbers,  a  reversal  of  trend  found  at 
low  8c«  This  is  due  to  the  rapid  increase  of  static  temperature 

on  the  cone  surface  with  increasing  6  .)  This  means  that  between 

c 

the  windward  surface  and  the  leeward  surface,  the  flow  goes  from 
subsonic  to  supersonic  conditions,  probably  through  a  sonic  line 
similar  to  the  flow  around  a  yawed  cylinder.  The  luminosity  of 
the  flow  seen  on  the  windward  surface  of  these  bodies  is  due  to 
slight  contaminants  in  the  flow. 


Figures  9  to  14  show  the  reduced  pressure  data  in  terms 

of  the  pressure  coefficient  c  ■  p/-^  pV2.  Here  the  dynamic  pres- 
1  9  P  2 

*ure,  —  pVz ,  at  the  nose  position  is  used.  The  results  are  shown 
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conpared  to  the  tangent-cone  and  Bevtonian  theorise  in  which  the 
corrections  outlined  ia  Appendix  C  are  carried  out.  The  value  of 
c  for  the  two  theories  for  a  50®  half-angle  cone  in  a  parallel 
flow  is  given  by  that  value  at  which  the  theories  cross;  *he  s»0 
ordinate.  The  results  are  seen  to  agree  within  a  few  percent  with 
the  tangent-con©  theory  which  itself  is  roughly  5%  above  Bevtonian 
theory.  Exceptions  are  seen  in  Runs  2G9  and  210  it  which  the  first 
and  third  gauges  give  low  values.  It  is  thought  that  this  is 
caused  by  non-linear  calibration  curves  (see  section  2.3).  The 
results  for  measurements  from  0  to  10  msecs  after  the  peak  are 
presented  in  these  figures. 

Fig.  15  show  the  heat  transfer  rates  at  peak  conditions 
plotted  against  local  Reynolds  nusber  for  the  sero  incidence  cases. 
Three  different  heat  transfer  predictions  are  compared  to  these 
results.  These  are  the  Eckert  reference  enthalpy  laminar  theory, 
the  Sommer  and  Short  reference  enthalpy  and  Spalding-Chi  turbu¬ 
lent  theories.  The  formulae  used  in  these  predictions  are  listed 
in  Appendix  C.  It  is  seen  that  all  the  theories  agree  reasonably 
well  with  the  experiments.  The  trend  however  is  best  predicted 
by  the  laminar  theory.  Alao  the  turbulent  theories  do  not  compare 
well  with  each  other.  The  flow  conditions  on  the  surface  of  the 
cone  as  shown  in  the  read  out  of  the  program  at  the  end  of  Appen- 
&•/  C  are  M«1.45,  T-»2280°K  and  Reynolds  number  based  on  the  gene¬ 
rator  length,  1.15x10s.  The  Reynolds  number  is  low  essentially 
because  there  is  a  large  increase  in  temperature  of  the  flow  over 
the  free-stream  condition  with  only  a  small  change  in  density. 

The  temperature  is  high  enough  that  the  gas  at  the  edge  of  the 
boundary  layer  is  in  a  vibrationally  excited  state.  Furtherac-re 
these  conditions  are  well  outside  the  range  of  conditions  that 
these  engineering  theories  were  correlated.  Even  though  the  turbu¬ 
lent  theories  are  as  close  to  the  experiment  as  the  laminar  theories, 
since  the  Reynolds  number  is  lew  and  the  bread  is  better  predicted 
b7  the  Sckert  theory,  the  boundary  layer  is  more  likely  to  be  in 
a  laminar  state.  Further  discussion  of  this  topic  is  jiven  in 
•action  h.3  involved  with  tests  on  the  blunt  cone. 


Figs.  1 6  and  17  show  the  effect  of  incidence  on  the 
heat  transfer  rate.  The  agreement  with  laminar  theory  is  fair  on 
the  "leeward"  surface  hut  poor  on  the  vindvard  surface.  However 
as  mentioned  earlier,  the  extrapolations  of  cone  tables  to  these 
conditions  are  particularly  severe  and  nakes  meaningless  any 
attempt  at  making  conclusions. 

4.2  Pointed  rough  biconic  model  (B) 

The  heat  transfer  distributions  on  the  rough  model  are 
shovn  plotted  in  Figs.  18-20.  It  is  seen  that  there  is  c.  consi¬ 
derable  scatter  of  data  along  the  model.  The  scattered  heat  trans 
fer  rates  and  pressure  measurements  (shovn  in  Fig.  21)  obtained 
vith  this  model  indicate  that  there  i«  appreciable  viscous-invis- 
cid  interaction  of  the  flov  over  the  sensors,  vhich  are  buried 
belov  the  tops  of  the  roughness  elements.  The  lov  values,  for 
instance,  could  occur  if  the  sensors  are  in  a  separation  region. 
Excluding  these  lov  values,  the  results  show  a  definite  trend 
tovards  the  turbulent  boundary  layer  theories.  This  is  particu¬ 
larly  seen  in  the  test  at  M»15  (Fig.  20)  in  vhich  the  laminar 
and  turbulent  theories  are  more  separated. 


Several  interesting  features  can  be  seen  from  the 
schlieren  pictures  of  Fig.  8,  Runs  211  to  214.  In  all  cases,  the 
bov  shock  remains  attached  to  the  nose.  Shock  vaves  emanating 
from  the  roughness  elements  can  be  seen  in  the  shock  region  of 
Run  2 1 4 .  This  M«15  case  gave  higher  density  levels  than  the  other 
cases,  alloving  the  gradients  to  be  seen  more  clearly. 

4.3  Blunt  smooth  biconic  model  (C) 

Fig.  22  shovs  the  reduced  pressure  data  in  terms  of  the 
pressure  coefficient  e^,  for  the  peak  values  of  each  of  Run  284 
to  Run  287.  Shovn  plotted  are  the  theories  of  Belotserkovskii 
(Ref.  11)  for  the  spherical  nose,  and  the  tangent-cone  theory  of 
Appendix  C  for  the  cone  surface.  A  dotted  line  shovs  the  probable 
pressure  variation  after  the  overexpansion  region  aft  of  the  nose 
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It  it  seen  that  the  pressure  recovers  from  the  overexpansion 
near  gauge  number  4.  The  data  for  the  zero  incidence  case  and 
the  leevard  surface  agree  very  veil  vith  the  predictions  shovn. 
However,  for  the  vindvard  case,  the  measurements  are  somewhat 
lover  than  the  tangent-cone  theory.  This  could  be  caused  by  the 
high  cross-flow  on  this  surface,  enhanced  by  the  subsonic  flow- 
field. 

Fig.  23  shows  the  reduced  heat  transfer  data,  non- 
dimensionaliaed  by  the  measured  stagnation  point  heat  transfer 
predictions  for  the  zero-incidence  cases.  Run  281*  and  287.  The 
results  are  compared  vith  the  laminar  similarity  theory  of  Lees 
(Ref.  3)  uncorrected  for  conicity  effects  and  also  the  Eckert 
reference  enthalpy  method.  In  these  cases,  the  experiments  agree 
well  vith  these  laminar  theories.  Also  the  stagnation  point  heat 
transfers  for  these  cases  agree  well  with  the  theory  of  Fay  and 
Riddell  (Ref.  2).  This  lends  support  to  the  supposition  that  the 
data  on  both  the  blunt  and  pointed  cones  are  in  fact  laminar. 

The  heat  transfer  data  for  the  yawed  cone  is  shovn  in 
Fig.  24  compared  vith  a  line  drawn  through  the  unyaved  cone  data 
of  Fig.  23*  It  is  seen  that  the  heat  transfer  rates  for  the  wind¬ 
ward  side  are  much  closer  to  the  zero  angle  of  attack  results 
than  those  measured  on  the  leevard  side.  This  is  a  similar  trend 
as  seen  in  the  pressure  distributions  of  Fig.  22. 

The  schlieren  photographs  in  Fig.  8  for  Runs  284  to  287 
show  similar  effects  to  that  seen  in  the  photographs  for  the 
pointed  body  except  in  the  now  detached  shock  at  the  nose.  Simi¬ 
larly  the  shock  shapes  suggest  that  supersonic  flow  is  found  on 
all  surfaces  except  the  vindvard  surface. 
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1*.i*  Hemisphere  models  (D.E) 

F.igure  25  shows  the  distribution  of  pressure  non-Aimen- 
sionalised  by  an  estimated  stagnation  pressure  for  all  fire  tests 
on  smooth  and  rough  hemispheres  and  compared  with  the  theory  of 
Belot3erkovskii  (Ref.  11)  corrected  for  eonicity.  It  is  seen  that 
agreement  is  excellent.  The  Pitot  pressures  have  been  plotted  as 
the  hemisphere*  stagnation  point  value.  It  is  seen  that  these  pres¬ 
sures  are  about  53?  higher  than  suggested  by  the  hemisphere  data. 
These  high  readings  may  be  explained  by  the  Pitot  probe  sensing 
the  pressure  in  the  ‘’overshoot"  region  of  the  nozzle  turbulent 
boundary  layer  region  ac  measured  in  Longshot  previously  in 
Ref.  12.  This  can  explain  th ^  high  degree  of  scatter  obtained  in 
the  measurements  of  Pitot  pressure. 

The  heat  transfer  measurements  for  the  smooth  hemisphere 
have  been  compared  with  the  similarity  theory  of  Lees  (Ref,  3) 
corrected  for  eonicity  and  with  a  rough  correction  to  align  the 
theory  with  the  pressure  curve  of  Belotserkovskii  (Ref.  11)  instead 
of  Bewtonian  theory  (Fig.  26).  The  measurements  were  non-dimen- 
sionalised  with  respect  to  the  theoretical  value  of  heat  transfer 
rate  of  Fay  and  Riddell  (Ref.  2).  Within  30  degrees  of  the  stag¬ 
nation  point,  the  measurements  are  seen  to  be  considerably  higher 
than  theory.  These  results  are  in  disagreement  with  earlier  mea¬ 
surements  on  a  wooden  hemisphere  model  at  M*15  reported  in  Ref.  8 
and  with  the  measurements  on  the  stagnation  point  of  model  C. 

A  possible  reason  why  the  wooden  model  results  gave  a  different 
result  than  in  these  tests  is  that  the  calorimeter  gauge  was  not 
contoured  to  the  surftce.  It  was  assumed  at  the  time  that  for 
these  large  model  radii,  the  disc  contouring  becomes  unimportant. 
However  this  discrepancy  io  puzzling.  The  differences  are  too 
great  and  are  too  coincidental  to  be  explained  by  uncalibrated 
gauges.  These  gauges  also  appeared  not  to  be  damaged  or  eroded 
during  a  shot.  Possible  suggestions  are  tunnel  contaminants  or 
free-stream  turbulence,  although  this  does  not  explain  why  such 
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good  results  were  obtained  with  the  0.75  in  radius  nose  of  model  C, 
After  30w,  then  the  results  agree  very  well  with  the  similarity 
theory  of  Leee,  thus  indicating  that  the  flow  is  laminar. 


Fig.  27  shows  the  results  of  the  rough  model  displayed 
in  a  similar  fashion  to  the  smooth  model  data.  Again  very  high 
stagnation  point  heating  rates  are  obtained.  High  heating  rates 
are  seen  also  generally  around  the  cone*  which  could  be  ascribed 
to  either  turbulent  flow  or  laminar  flov  undergoing  mixing  by 
the  roughness  without  going  turbulent.  The  highest  relative  heat 
transfer  rates  are  found  in  the  M**15  case.  This  is  the  mo?t  likely 
case  that  turbulent  flow  could  be  promoted. 

The  shock  staed-off  distances  from  the  nose  of  models 
C,  D  and  E,  as  measured  from  the  schlieren  pictures  of  Fig.  8 
are  shown  plotted  in  Fig.  28  and  compared  with  the  theory  of 
Van  Dyke  as  presented  in  Ref.  19.  The  comparison  is  ahown  to  be 
very  good. 
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liJIM  tiflw  i 


15 


5.  C0HCLU8I0I3 


Tests  have  been  carried  out  on  3  biconic  and  2  hemi- 
sphere  nodels  in  an  M«15  and  20  flow  in  the  VKI  Longshot.  The 
data  has  been  compared,  preliminarily,  vith  simplified  theories 
to  assess  their  quality.  The  following  conclusions  can  be  made 
from  the  comparisons. 

1.  Pressure  data  un  50°  aones  and  hemispheres  agree 
very  veil  vith  tangent-cone  and  the  Belotserkovskii  theory  res¬ 
pectively.  An  exception  vas  found  on  the  vindvard  side  of  the  50° 
blunt  cone  when  lover  pressures  than  expected  vere  found.  This 
wits  probably  caused  by  cross-flow  enhanced  by  locally  subsonic 
flov  behaviour,  or  due  to  the  severe  extrapolation  used  in  the 
theory. 


2.  High  quality  and  spatially  smooth  heat  transfer 
data  was  obtained  on  all  models  except  the  very  rough  biconic 
model.  There  is  some  doubt  as  to  whether  the  cone  data  is  laminar 
or  turbulent  due  to  the  similarity  in  predictions  of  laminar  and 
turbulent  theories  at  the  Longshot  flow  conditions  at  M«20.  The 
evidence  weighs  more  heavily  vith  the  presence  of  laminar  flov 
due  mainly  to  the  better  prediction  in  trend  of  the  experimental 
results  by  the  Eckert  reference  enthalpy  method.  Heat  transfer 
measurements  on  the  rough  models  give  good  indications  that 
transition  to  turbulent  flows  has  occurred,  especially  in  the 
K*15  cases. 

3.  Stagnation  point  beat  transfer  on  the  blunt  cone 
agreed  very  veil  with  Fay  and  Riddell,  but  some  discrepancies 
have  been  found  vith  the  7  is  hemisphere  results. 

fc.  Sfc-  c'y  shapes  on  the  conical  models  and  shock  stand¬ 
off  distances  from  the  spherical  caps  of  all  the  blunt-nosed 
models  agree  very  veil  vith  predictions. 


A  general  conclusion  in  that  useful  laminar  and  turbu¬ 
lent  boundary  layer  heat  transfer,  pressure  and  flov  visualisa¬ 
tion  data  has  been  achieved  on  nose  shapes  in  a  flow  simulating 
aerodynamic  re-entry  conditions. 
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APPEKDIX  A 


A  Short  Description  of  the  Longshot 
Data  Reduction  Programme 

A  computerised  data  reduction  programme  has  >een  devised 
by  Richards,  Culotta  and  Enkenhus  and  is  in  general  use  in  Longshot. 
The  details  vill  be  written  up  in  a  future  report,,  It  is  appro¬ 
priate  to  give  a  brief  description  of  the  method  in  this  report 
for  understanding  the  fesults, 

In  attaining  very  advanced  test  conditions,  which  simu¬ 
late  re-entry  flight,  the  V.K.I,  Longshot  free-piston  wind  tunnel 
displays  some  unusual  features.  During  a  test,  dense  real  nitrogen 
is  trapped  in,  and  allowed  to  escape  from,  a  constant  volume 
reservoir  into  a  conical  nossle  to  provide  a  high  Reynolds  number, 
high  Mach  number  test  flow.  It  has  been  shown  theoretically  and 
experimentally  (Ref,  13)  that  the  reservoir  and  test  section  para¬ 
meters  measured  in  Longshot  can  be  described  accurately  by  an 
exponential  decay  function,  f,  of  the  form  t 

f  »  exp  (ao  +  ajt  +  ajt2)  A.l 

where  t  is  the  time  in  milliseconds  from  the  peak  value  of  f. 

All  the  measured  reservoir  and  test  section  pressures 
and  heat  transfer  rates*  are  curve  fitted  to  this  formula  using 
a  least-squares  method  (Fig,  A.l).  The  differences  between  the 
measured  and  re-generated  values  using  equation  A.l  are  found 
to  be  within  reading  error. 

The  supply  temperature  is  measured  with  a  thermocouple 
with  a  long  time  response*.  Use  is  made  of  a  formula  which  connects 
the  temperature  explicitly  with  supply  pressure  (see  Ref,  13) 


Actual  traces  of  reservoir  pressure  and  temperature  and  Pitot 
and  wall  pressures  are  shown  in  Fig,  5, 
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as  it  decays  fro®  a  constant  volume  reservoir  to  calculate  the 
temperature  variation  backwards  in  time  from  a  time  when  the 
thermocouple  is  assessed  to  be  giving  the  true  temperature 
(Fig,  A, 2),  Conduction  losses  from  the  wire  and  radiation  losses 
from  the  gas  are  taken  into  account  to  assess  the  true  supply 
temperature.  This  difficult  temperature  measurement  will  be  the 
object  of  a  further  future  report  by  Richards  and  Enkenhus, 

The  test  gas  expands  from  the  reservoir  until  it  becomes 
a  perfect  gas  in  the  test  section,  The  density  of  the  gas  during 
the  early  stages  of  expansion  is  high  enough  that  freezing  of  the 
internal  energy  modes  does  not  occur  and  thus  the  expansion  can 
be  considered  as  being  in  equilibrium,  A  concept  of  equivalent 
perfect  supply  conditions,  described  in  Ref,  }  is  used  to  faci¬ 
litate  the  calculation  of  the  test  flow  parameters. 

Supplying  the  basic  recorded  oscilloscope  or  oscillo¬ 
graph  deflections  and  instrumentation  calibration  data  from 
Longshot's  fast  response  instrumentation,  the  programme  can  cal¬ 
culate  supply  pressure  and  temperature,  Pitot  pressure  and  sur¬ 
face  heat  transfer  and  pressure  on  models,  Xf  the  two  supply 
conditions  end  a  test  section  Pitot  pressure  is  supplied,  the 
programme  calculates  the  free  stream  conditions  including  Mach 
number,  Reynolds  number  and  dynamic  pressure  (as  described  in 
Ref,  14) ,  the  stagnation  point  heat  transfer  (Ref,  1  and  as  in 
Appendix  D)  and  the  test  gas  condensation  temperature  (Ref,  15), 

A  typical  set  of  data  obtained  in  this  present  series 
of  tests  is  shown  cn  the  next  few  pages. 

The  first  part  of  tbe  programme  for  each  measurement 
gives  the  calibration  constants  and  the  trace  deflections,  Y 
divisions  (cm  or  mm)  at  equal  time  intervals  (as  indicated  by  T) 
from  an  easily  identifiable  position  on  tbe  trace,  Y  RAW  gives 
the  raw  data,  Y  CORR  the  smoothed  data.  Except  in  the  case  of 
the  reservoir  temperature,  the  difference  between  Y  RAW  and  Y  CORR 
is  seen  to  be  very  small.  Comparison  of  Y  RAW  and  Y  CORR  under 


reservoir  temperature  illustrates  the  extrapolation  method,  used* 
The  values  agree  at  the  later  tines  in  the  rune,  when  the  thermo- 
couple  has  reached  its  equilibrium  temperature,  ’-at  deviate 
vastly  in  the  early  times  shoving  the  slov  response  of  the  thermo 
couple*  This  difficult  measurement  is  the  least  accurate  of  those 
made  in  the  Longshot ,  but  is  considered  to  be  vithin  an  error  of 
±5*. 

The  second  part  of  the  programme  ("quantities  measured 
in  Longshot")  summarises  all  the  measurements  at  times  from  the 
peak  measured  conditions  using  the  appropriate  curve  fit  of  the 
actual  data  to  Equation  A.l,  If  a  Pitot  pressure,  p.  ,  is  given 
then  for  the  vail  pressures  an  extra  rov,  YP  is  given  vhich  gives 
the  pressure  non-dimens ionalized  to  the  dynamic  pressure,  1/2  pV2 
vhich  is  itself  approximated  using  the  formula  : 

1/2  pV2  -  0.5fc3  Pt 

z  2 

This  formula  is  true  for  a  perfect  gas  at  hypersonic  speeds,  but 
the  Pitot  pressure  measured  in  Longshot  is  different  from  a  per¬ 
fect  gas  value  by  about  ljt  (the  actual  correction  is  given  in 
Fig*  of  Ref,  }  due  to  the  gas  being  vibrationally  excited 
behind  the  normal  shock  vave  formed  ahead  of  the  Pitot  probe. 

Thus  Y?  is  approximately  IJt  lover  than  the  true  value.  This  is 
veil  vithin  the  limits  of  the  measurements  and  since  it  is  used 
only  as  a  rough  guide,  it  has  not  been  thought  necessary  to  make 
the  correction. 

The  third  part  of  the  programme  ("Longshot  test  condi¬ 
tions")  reads  out  at  the  time  intervals  from  the  peak  the  follov- 
ing  variables  :  measured  reservoir  pressure,  PO  (psi);  measured 
reservoir  temperature,  TO  (°K);  measured  Pitot  pressure,  PITOT 
(psi);  calculated  Mach  number;  equivaleut  perfect  pressure,  POP 
(psi);  equivalent  perfect  temperature,  TOP  (°K);  Reynolds  number 
per  ft;  local  freestream  pressure,  F  (psi);  freeetream  tempera¬ 
ture,  T  (°K);  freestream  density,  RHO  (slugs/ft3);  stream  velo¬ 
city,  V(ft/see);  dynamic  pressure,  QD  (lb/ft2);  stagnation  point 


-  A  !»  - 

heating  on  a  7  in  diameter  spherical  surface,  Q  (Btu/ft2sec) j 
true  stagnation  temperature,  TT2R  (°K)i  and  the  teaperature  at 
which  condensation  would  occur  at  that  free-stream  pressure  and 
expansion  rate,  COHDEHSATIOH  (°K),  These  parameters  adequately 
define  all  the  parameters  necessary  for  application  to  predictive 
procedures.  Care  should  be  taken  in  that  the  accuracy  of  the 
values  printed  out  should  not  be  inferred  from  the  six  signifi¬ 
cant  figures  shown.  The  accuracy  is  controlled  by  the  accuracy 
of  the  measurements  .inferred  in  section  2. 
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LONGSHOT  TEST  CONDITIONS,  RIJN  NO. 286  AT  7.5,  IN.  .FROM  NOZZLE  EXIT, 
AND  4.6,  FROM  CENTRELINE 


T(MS) 

MACH  NO 
P(PSI) 

QD(LB/FT**2) 

POCPSf ) 

POP(PSI) 

T(K) 

Q(BTU) 

TD(K) 

TOP(K) 

RHO 

TT2R(K) 

PlTOT(PSI ) 
RE/FT 

V (FT/SEC) 
CONDENSATION 

0.0,00 
19.046 
0,1.494 10E-01 
0.53918QE  03 

0.5,9308 2 E  05 
0.7.01402E  05 
0.3.87404E  02 
0.6.79880E  02 

0. Z41827E  04 
0.312127E  04 

0. 173726E-04 
0.2S7115E  04 

0.776414E  01 

0. 291117E  07 
O.8.3054OE  04 
0.356669E  02 

2,000 
19.  7.86 
0.1.40848E-01 
0.5.55816E  03 

0. 5.35447E  05 
0.625396E  05 
0.3.76368E  02 
0.621384E  02 

0 . 233888E  04 

0. 298450E  04 
0.168572E-04 

0. 252752E  04 

0, 7.21249E  01 
0.286180E  07 
0.812057E  04 
0.354772E  02 

4.0,00 
19.053 
0.132779E-01 
0.5,169 8 7E  03 

0.48 705 IE  05 
0.5.62846E  05 
0.366010E  02 

0. 5.7099 5E  02 

0.226723E  04 

0. 286415E  04 

0 . 163412E-04 

0. 241631E  04 

0.S71268E  01 
0.281271E  07 
0.795448E  04 
0.352868E  02 

e.ooo 

19.5.48 

0.125120E-01 

0. 4.81985E  03 

0.446367E  05 
0.511142E  05 
0.356269E  02 
0.5.27424E  02 

0.220278E  04 
0.2758S8E  94 
0.158201E-04 
0.232612E  04 

0.6.25931E  01 
0.276314E  07 
0.780597E  04 
0.350941E  02 

8.0100 
19.4.72 
0.117814E-01 
0.430297E  03 

0.412161E  05 
0.468438E  05 
0.3.46957E  02 

0. 4.89543E  02 

0. 2145Q2E  04 

0. 266584E  04 
0.152957E-04 
0.225048E  04 

0.584758E  01 
0.271408E  07 
0.7.67324E  04 
0.3.48982E  02 

9.R99 
19.427 
0. 1.10  7  99  E- 01 
0.4.21530E  03 

0.383440E  05 
0.433561E  05 
0.337869E  02 
0.456433E  02 

0. 209354E  04 
0.258418E  04 
0.147719E-04 
0.218561E  04 

0.547325E  01 
0.266766E  07 
0.755459E  04 
0.346977E  02 

12.000 
19.408 
Q.L04120E-01 
0.3.95562E  03 

0.359409E  05 
0.404735E  05 
0.3.29141E  02 
0.427462E  02 

0.204792E  04 
0. 251266E  04 
0. 142496E-04 
0. 212958E  04 

0.513256c  01 
0.262216E  07 
0. 744922E  04 
0.344944E  02 

14.Q00 
19.415 
0.9.77699E-02 
0.3.71521E  03 

0.3.39418E  05 
0.3,8 10 18 E  05 
0.320768E  02 
0.402118E  02 

0.200786E  04 
0. 245053E  04 
0.137297E-04 
0. 208123E  04 

0. 482216E  01 

0 . 257739E  07 
0.735657E  04 
0.342885E  02 

16.Q00 
19.448 
0.9.17379E-02 
0.3A9775E  03 

0.3.22954E  05 
0.3,6 169 IE  05 
0.3.12737E  02 
0.3.79965E  02 

0.197306E  04 
0. 239712E  04 
0.132135E-04 
fj .  203977E  04 

0.4.53910E  01 
0.253328E  07 
0.727612E  04 
0.3408«1E  02 

17.939 
19.006 
0. 8.60166E-02 
0.339918E  03 

0. 3.096 03E  05 
0.3.46176E  05 
0.3.95039E  02 
0.3.60637E  02 

0.194327E  04 
0.235189E  04 
0.127020E-04 
0. 209466E  04 

0.4.28072E  01 
0.248973E  07 
0.720743E  04 
0.338695E  02 

H.T,.  DECAY  COEFF...LAM.  .  -0.4748E-01  0.4800E-03 

TURB.  -0.4995E-01  0.4794E-03  STAG  -9.4.5606-01  0.5.B07E-03 


APPESDIX  B 


Data  Reduction  of  Calorimetric  Heat  Transfer  Data 

A a  in  the  previous  appendix,  this  data  reduction  pro¬ 
gramme  has  been  perfected  in  the  last  fev  months  this  time  by 
Culotta  and  Richards  and  will  again  be  the  subject  of  a  future 
report. 

Calorimeter  gauges  are  constructed  of  a  thin  disc  of 
copper  bonded  to  an  insulated  holder.  Thermocouples  vith  a  vire 
diameter  smaller  than  the  thickness  of  the  disc  are  bonded  (either 
welded  or  soldered)  to  the  back  face  of  the  copper  disc.  The 
thermocouple  senses  the  average  temperature  rise  of  the  disc  as 
it  heats  in  the  test  flov.  The  heat  transfer  rate,  q,  can  then 
be  determined  from  the  following  equation  : 


q  -  k 


dV 
c  dt 


B.l 


where  k  is  the  gauge  calibration  constant 
c 

dV/dt  is  the  slope  of  the  thermocouple  output  variation 
with  respect  to  time* 

The  calibration  constant,  kfi,  is  best  determined  expe- 
rimentally  using  a  rapidly  applied  and  calibrated  heat  source. 
However,  it  can  be  calculated  from  knowledge  of  the  calorimeter 
thickness  (b),  the  density  (p^)  and  specific  heat,  (c^)  of  the 
disc  material  and  the  thermocouple  response  parameter  dT/dV, 
using  the  following  equation  : 


k  *  b  p. 
c  d 


.  dT 
cd  dV 


B.2 


It  must  be  assumed  however  that  there  are  no  heat 
conduction  losses  to  the  holder  or  to  the  thermocouple,  the  disc 
thickness  is  not  changed  by  the  welding  or  soldering  of  the  thermo¬ 
couple,  the  temperature  is  uniform  over  the  volume  of  the  disc* 
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The  order  of  magnitude  of  the  tiae  constant,  t,  for 
the  calorimeter  may  be  assessed  from  : 

Pde.b; 

b.3 

vhere  k  is  the  theraal  conductivity. 

For  the  testing  of  a  model  at  rest  in  a  steady  flov, 
equation  B.l,  is  simple  to  apply  for  calculating  the  heat  trans¬ 
fer  rate  from  the  raw  thermocouple  output  data.  The  heat  trans¬ 
fer  rate  would  be  expected  to  be  constant  and  thus  the  output 
will  vary  linearly  with  time.  It  is  a  simple  matter  to  make  a 
linear  least  error  fit  to  the  data  points  to  achieve  the  measure¬ 
ment.  The  method  depends  upon  an  assumed  knowledge  of  the  vari¬ 
ation  (but  not  the  level)  of  the  heat  transfer  with  time,  in  this 
case  it  remains  constant. 

For  varying  flov  conditions  alternative  methods  must 
be  found.  More  usual  methods  are  to  make  a  polynomial  curve  fit 
to  the  raw  data  and  then  to  differentiate  the  formula  thus  ob¬ 
tained.  Often  poor  results  are  obtained.  If  a  large  order  poly¬ 
nomial  is  chosen  although  an  accurate  curve  fit  to  the  data  is 
achieved,  the  derivative  can  fluctuate  wildly  about  its  true 
variation  due  to  very  small  errors  in  reading  the  data.  If  a  small 
order  polynomial  is  chosen  to  eliminate  these  effects,  there  is 
a  risk  of  losing  important  detail.  In  general,  differentiating 
a  polynomial  curve  fit  to  data  is  usually  unsatisfactory. 


A  new  data  reduction  technique  has  been  devised  which 
uses  an  extra  piece  of  information  to  allow  the  reverse,  but  more 
accurate  process  of  the  integration  of  equation  B.l,  The  extra 
information  is  the  knowledge  of  the  variation  (but  again  not  the 
level)  of  the  heat  transfer  rate,  given  in  the  form  of  Equation 
A.l,  calculated  from  the  calibrated  test  section  from  the  Longshot 
data  reduction  programme  (Appendix  A),  Three  engineering  calcu¬ 
lations  of  the  heat  transfer  rate  (i.e.  the  laminar  and  turbu¬ 
lent  heat  transfer  rate  at  a  fixed  position  on  a  flat  plate  and 


. . 


the  stagnation  point  heating  rate  on  a  sphere),  are  made  in  this 
latter  programme,  and  the  coefficients  a*  and  a*  of  equation  B,1 
are  determined,  These  are  read  out  at  the  very  end  of  the  prograaae 
as  shown  in  Appendix  A.  Trial  values  of  the  heat  transfer  at  tero 
tiae  (4^,0  *  ft«0  *  **°)  're  *ade  and  equation  B,1  is  integrated 
vith  this  coaplete  inforaation  for  q.  The  calculated  transient 
teaperature  variation  thus  obtained  is  compared  vith  the  measured 
thermocouple  output.  The  value  of  q^tJgQj  that  provides  the  best 
comparison  of  the  temperature  variations  is  considered  the  correct 
value.  The  remainder  of  the  variation  of  q  vith  tiae  is  then  cal- 
culated  and  read  out.  The  numerical  integration  is  carried  out 
using  a  Runge-Kutta-Gill-Simpson  routine  and  the  search  proce¬ 
dure  used  converges  very  rapidly. 

The  above  technique  is  analogous  to  th?  ssihod  described 
above  for  dealing  vith  constant  heat  transfer  data*  To  suuaarise, 
in  both  methods  an  assumption,  in  made  about  the  variation  (but 
not  the  level)  of  the  required  parameter  ,  in  this  case  the  heat 
transfer  rate.  There  is  no  reason  vhy  this  method  cannot  be  ex¬ 
tended  to  data  achieved  from  more  complicated  unsteady  test 
section  conditions  than  in  Longshot  vhen  an  appropriate  equation, 
analogous  to  equation  A.l,  must  be  used  t-o  describe  the  heat 
transfer  variation. 

To  illustrate  the  technique  a  typical  output  from  the 

computer  is  shovn  later.  A.l  and  A. 2  are  the  coefficients  defining 

the  slope  and  curvature  data  as  read  from  the  last  line  in  the 

Longshot  data  reduction  programme  for  the  appropriate  case  (i.e. 

stagnation,  laminar,  turbulent  heat  transfer  case).  K1  and  K2  are 

the  oscilloscope  calibration  in  mV/cm  or  mV/mm  and  the  gauge 
B  Th  0  mV 

calibration  in  — *"■*■  /  — *  .  QO  is  the  value  of  heat  transfer 
ft2sec  8ec 

rate  in  B.Th.U/ft2sec  at  zero  time  at  each  iteration  step  until 
convergence  is  obtained.  (It  is  seen  that  convergence  is  reached 
almost  immediately).  Q(T1)  is  the  final  value  of  heat  transfer 
rate  at  zero  time  (peak  value).  T  is  the  time  from  peak  condi¬ 
tions  in  msecs,  Y  EXP  is  the  measured  oscilloscope  deflection 
(cm  or  mm).  Y  CALC  is  the  equivalent  calculated  value  vith  the 


ttlilll _ _ _ _ 
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best  fit  to  the  measured  data  a*  outlined  above,  ERR  is  the 
difference  between  Y  EXP  and  Y  CALC.  (This  value  iB  generally 
within  reading  errors  as  can  be  seen  in  the  example).  Q  is  "whe 
heat  transfer  rate  <3.Th.U./ftzsec)  at  time  T.  For  Q(Tl)  then 
since  the  calculation  starts  at  T  ■  0,  then  ERR  *  0, 


DATA  REDUCTION  OF  CALORIMETER  HEAT  TRANSFER  GAUGES,  RUN  286,  GAUGE  0 
Al*-0.4.9980E-01  A2«  0.48050E-03  Kl*  0.10100E  00  K2a  0.73200E  00 

Q0»  0.149792E  03 

Q0-  0.149792E  03 

CONVERGENCE 

Q(T1>  =»  0.14979E  03 


T 

Y  EXP 

Y  CALC 

ERR 

Q 

2.000 

0.52000E 

01 

0.52586E 

01 

0.58643E-Q1 

0.13580E 

03 

4.0.00 

0.9.100QE 

01 

0.87635E 

01 

-0. 33643E  00 

0.12359E 

03 

6.000 

0.12000E 

02 

0.11959E 

02 

-0.40557E-01 

0.11291E 

03 

8.000 

0.14900E 

02 

0.14884E 

02 

-0.15230E-01 

0.10355E 

03 

10.0.00 

0.17500E 

02 

C.17572E 

02 

0.72757E-01 

0.95345E 

02 

12.000 

0.20000E 

02 

0.20052E 

02 

0.52177E-01 

0.88118E 

02 

r 


\ 


1 


1© 


I 
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API-EftDIX  C 

i 

Prediction  of  Beat  Transfer  Hate 
on  a  Sharp  50°  Half-Angle  Gone 

i 

A  computer  programme  vas  written  by  Richards  and 
Slaehten  to  estimate  the  laminar  and  turbulent  flow  over  a  50° 
half-angle  cone  in  a  hypersonic  conical  flew.  The  equations  used 

t 

are  described  below, 

1.  Calculation  of  the  test  section  conditions 

The  initial  data  of  the  equivalent  perfect  supply 
conditions  (Tq)^^  an<*  ^D0^perf  and  the  Mach  number  at  an  arbi¬ 
trary,  but  knovnposition,  in  the  tunnel  are  supplied  as  basic  . 
input  data  to  the  programme.  These  values  are  obtained  from  the 
Longshot  data  reduction  programme* 


The  test  section  conditions  can  then  be  obtained  using 


the 

perfect  gao  relations  : 

• 

Tl 

-  T0(l  +  M2/5)“l 

(°K) 

C.l 

*1 

-  P0(l  +  M2/5)~7/2 

(lb/in2) 

C.2 

Pi 

*1 

"  rt7 

(slug/ft3) 

C,3 

*1 

*  /ySTi 

( ft/sec) 

C.b 

ui 

■ 

(ft/sec) 

C,5 

Ml 


93.52T  Tt3/*  a  .. 

■  -  x  in*3  f I, .  .1.^ ■ . ) 

Ti  +  102, 7  i0  'ftsec' 


C,6 


Rei 


(ft-1) 


i 


i 


i 


I 
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2*  Calculation  'of ‘inviscid  flov  over  a  cone  , 

The  cone  tables  .of  Jones  (Ref.  9  )  for  y-!*1*  vere  used 
to  determine  the  conditions  on  a  cone.  For  various  angles  of 

.  t 

attack,  the  tables  give1 the  parameters  : 


u  a 

c  c  c  ____ 
*1*  61*  Pi*  1 

i  ■  2 


PjUi2  sin2  6 


as  functions  of  M  and;0  .  Here  the  subscript  c  denotes  the  con- 

,  .  '  c  ,  •  ' 

ditions  on  the  cone  surface  and  subscript  s  the  conditions  at 
the  shock.  The  range  of  applicability  of  Ref.  is  1.5  <  M  <  20 

and  5°  <‘ 6  <  U0°.  This  means  that  for  the  model  under  test, 

c  .  , 

extrapolations  had  to  be  made  from  0  ’ *  U0°'to  50°. 

„  !  : 

i  The  tabulated  values  of  the  parameters  listed  above 

i  i  4  ,  i 

vere  fitted  to  a  quadratic  ;pover  series  vith  M  and  8c  as  the 
independent  variables.  The  power  series  was  : 

i  .  * 

'f  -  Al+A28c-*-A30|+(A^A5ec+A6e2)M+(A74A8ec+A9e2)M2  ,  C.8 

The  coefficients  Ay  ...  Ag  were  obtained  by  initially  fitting 
quadratics  to  three  sets  of  data  in  one  of  the'  independent 

i  * 

variables.  The  latter  coefficients  thus  obtained  vere  then  cross 
fitted  to  quadratics  in  the. other  independent  variable  to  obtain 

;  *  I 

the  final  values  of  the  coefficients,  k\  ...  A9,  The  technique 
is  in  effect  a  two-dimensional  interpolation  procedure. 

•  i 

!  i 

Three  subroutines  vere  made  each  for  the  case  of  a*0 

,  •  I 

and  at* ±10°  (called  SLE1,  8LE2  and  SLE3  respectively). 


.  The  pressure  on  the  cone  as  calculated  in  this  fashion 
gives  the  tangent  cone  value,  The  modified  Hevtonian  value  of 
pressure  is  given  by  piui,,2sin2(0+a) . 


WOTWWyW'Wwr*  . . . . . . . . . . . . . . . . 
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3*  Calculation  of  conicity  effects 

Two  primary  corrections  arc  necessary  to  account  for 
conicity.  One  is  based  on  one-dimensional  flow  considerations, 
the  other  on  two-dimensional  flow  considerations. 


Regarding  the  first  correction,  consider  a  hypersonic 
conical,  perfect  gas  flow.  The  energy  equation  for  this  flow  is 

2 

h0  ■  cpTj  +  ui/2  C.9 

where  ho  is  the  reservoir  enthalpy  and  c^  is  the  specific  heat 
at  constant  pressure. 


For  hypersonic  flow  is  very  small  compared  to 

h0.  and  thus  u*  becomes  nearly  constant.  Kost  of  the  internal 
energy  has  been  converted  to  kinetic  energy,  and  the  change  in 
Mach  number  comes  mainly  from  a  change  in  7; .  hence  e 
change  in  aj . 


The  continuity  equation  (pjUiA  ■  constant,  where  A  is 
a  stream  tube  area)  then  gives  an  expression  for  the  density 
variation,  p ,  at  any  streamwiae  position,  x.  in  a  conical  flow 
field  : 


.-l 


PI 


where  a  «  ~  ■  ( 1  +  £—) 
Ai  Lq 


C.10 


C.ll 


here  subscript  1  refers  to  conditions  at  a  reference  position. 

L  is  the  distance  from  the  source  point  to  this  position, 
n 


Assuming  an  isentropic  expansion  :  p/pj  *  o”Y 
Equation  of  state  :  T/Ti  ■  o 


C.12 


¥ 


Sound  speed  :  a/al  ■  a 


C  ,13 
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Mach  number  i  H/Mi  *  o 


I=i 

2 


c  *  lU 


Equations  C.10  to  C.l4  can  then  be  used  to  calculate  the  stream 
conditions  at  any  longitudinal  position  in  the  expanding  flow, 
knowing  the  conditions  at  another  position*  Refs*  and 
have  shovn  that  the  Longshot  flow  in  the  nozzle  is  truly  conical 
with  its  source  point  located  at  the  nozzle  throat. 

The  second  correction  factor  to  account  for  conic ity 
arises  since  the  flow  angle  changes  by  an  angle.  do  *  tan“l  z — 

Ii  tX 

at  a  distance  y,  from  the  centre-line  of  the  nozzle, and  x  dowi- 
stream  of  the  reference  plane. 

The  application  of  these  two  corrections  to  determine 
the  conditions  on  an  arbitrary  position  on  a  cone  are  carried  out 
using  the  interpolation  equations  (C.8)  from  the  cone  table  data 
using  the  following  information  for  the  basic  parameters  M  and 
6c  and  for  the  other  freestream  parameters.  The  freestream  con¬ 
ditions  (including  M)  are  taken  as  those  that  would  exist  at 
the  same  position  in  the  test  section  if  the  model  were  removed. 
The  cone  angle  to  be  used  to  apply  the  interpolation  equations 

is  taken  to  be  6  -La.  (For  a  model  at  zero  incidence,  referred 

c 

to  the  nozzle  centre-line,  the  correct  application  of  the  tables 
should  be  by  using  the  values  obtained  on  a  cone  of  half-angle, 

6  ,  at  an  angle  of  attack  do.  However  little  improvement  would 

C 

be  made  for  this  interpretation,  since  da  is  small  compared  to 

0  ).  The  method  has  been  shown  to  work  well  in  the  results  and 
c 

discussion  of  this  report.  The  success  is  due  to  the  closeness 
with  which  the  bow-shock  lies  to  the  surface  of  the  cone. 


4.  Heat  transfer  rate  predictions 
Laminar 

The  Eckert  reference  (intermediate)  enthalpy  method 
is  used  (Ref.  4).  The  reference  temperature,  T*  is  given  by  : 


T  ♦  0.5  (T  -  T  )  +  0.22  (T  -  T  ) 
c  w  c  r  c 


C.15 


The  equivalent  perfect  value  of  T  is  used  in  the  cal- 

r 

culation  and  hence  the  application  of  equation  C,15  is  equivalent 
to  applying  the  reference  enthalpy  method,  Ty  is  obtained  from  : 

T  ■  T  +  r  u2/2e  C.l6 

r  c  ,  c  p 

vhere  r  is  the  recovery  factor  having  a  value  of  P  °»s  and  Pr 

r 

is  the  Prandtl  number  taken  to  have  a  value  of  0.715, 

The  incompressible  skin  friction  formula  used  was  the 
BlaBius  equation  : 

c.  «  0.664//Re"  C.17 

X  8 

vhere  c^  is  the  local  skin  friction  and  Res  is  the  Reynolds 
number  at  station  s. 


The  heat  transfer  rate  can  then  be  calculated  from  : 


°\CF  T-  \  (Tr  '  V 


C.18 


vhere  is  the  density  at  the  edge  of  the  boundary  layer  calcu¬ 
lated  at  the  reference  temperature  T*,  u  is  the  velocity  at 
the  edge  of  the  boundary  layer,  kR  is  the  Reynolds  analogy  para¬ 
meter,  and  (Tr  -  Ty)  is  the  driving  temperature  for  the  equation, 
c^  is  the  Bkin  friction  at  reference  temperature  conditions. 


For  laminar  flov,  kR  is  normally  taken  an  Pr -2/3 
and  ths  skin  friction  coefficient,,  is  calculated  from  equation  17* 
In  this  equation,  the  Reynolds  number  is  given  by  : 


Re.  "  r-  +  -nr  Lb)  C.19 

8  l  8  n-1  y  f 

c 

k^  is  the  Mangier  transformation  factor,  used  to  transform  a  skin 
friction  lav  from  a  tvo-dimensional  to  an  axisymmetric  case  (tbe 


i 
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flow  over  an  axisymmetric  body  has  a  thinning  action  on  the 
boundary  layer) »  k^  has  a  value  of  3  for  laminar  flow* 

The  term  in  the  brackets  is  an  accumulative  Reynolds 

number,  used  in  this  varying  flow  field.  (Re  )  is  the  aceu- 

8  n-1 

mulative  Reynolds  number  at  the  last  computation  step.  As  is 

the  step  size,  p*  and  y*  are  the  density  and  viscosity  of  the 
c  c  . 

flow  on  the  cons  calculated  at  the  reference  temperature, 
is  given  by  equation  C.6. 


Turbulent 

Tvo  theories  were  used  :  the  Sommer  and  Short  refe* 
rence  enthalpy  method  (Ref.  IT  )  and  the  Spalding-Chi  method 
(Ref.  5). 


The  Sommer  and  Short  method  vas  applied  in  a  very  si- 
milar  fashion  to  the  laminar  Eckert  reference  method.  However 
the  Sommer  and  Short  reference  temperature  is  given  by  : 


■  T  (1  +  0.035  M2  +  0.45  (T  /T  -  l)) 
c  w  c 


C  .20 


The  turbulent  recovery  temperature  is  given  by  : 


T  «  T  +  (Pr)1/3  u2/2c 
r  c  c  p 


C  .21 


The  incompressible  skin  friction  lav  used  vas  the 
Karmr.n-Schoenherr  implicit  equation  (Ref.  )  : 


,  _  0.242 

:f  log( Re  *  cf' 


0.557  cf> 

and  c.  *  — - - - 

0.557  +  2(cf')1/2 


C  .22 


C  ,23 


c^'  is  the  average  skin  friction  coefficient  and  is  the  local 
skir  friction  coefficient.  The  above  simple  equation  is  a  best 
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fit  to  the  universal  skin  friction  curve. 


There  is  such  uncertainty  shout  the  value  of  the 
Reynolds  analogy  factor,  kR  (see  Ref,  )*  An  average  value  of 


kD  *  1,08  vas  chosen  for  this  study.  In  the  turbulent  version 

X\ 


of  equation  C.19,  the  Mangier  transformation  factor,  kg,  has  a 
value  of  2  (Ref.  4). 


The  Spalding-Chi  method  (Ref.  5)  vas  applied  by  ob¬ 


taining  the  localskin  friction,  cf,  from  the  tabulated  data 


given  in  Ref.  5  using  input  values  of  Re  (turbulent  version  of 


equation  C.19),  M  and  T*  *  T  /T  ,  A  quadratic  series  vas  used 

0  VC 


to  calculate  C*(Re  ), 
r  s 


The  heat  transfer  equation  used  vas  ; 


(,,r  -  Tv> 


am 
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HEAT  TRANSFER  RATE  ON  CONE,  RUN  2C5,  AT  0.0.  MSECS  FROM  PEAK. 

ALPHA-  0.0  DEGREES.  TIP  AT  0.00 IN  ABOVE  CENTRE  LINE 

EQUIVALENT  PERFECT  CONDITIONS 
TO  -  0.350483E  04  DEG  K 

PO  «•  0.6A5708E  05  PSI 

M  -  0.20D050E  02 

THETC  -  0.5<0000QE  02  DEG 

TV/  -  0.23L.Q02E  03DEG  K 

PR  -  0.7.15000E  00 

FREE  STREAM  CONDITIONS  (ID  AT  TIP  OF  CONE 

T  -  0* 195462E  02  DEG  K 

P  -  0.1A2700E-01  PS! 

RHO  -  0, 162544E-04  SLUGS/FT**3 

U  -  0. 8416G7E  04  FT/SEC 

MU  -  0.5O8279E-07  LB*SEC/FT**2 

RE  -  0.2JS9140E  07 

CONE  CONDITIONS  AT  TIP  (EQUIVALENT  PERFECT) 

TE  «  0.2278KE  04  DEG  K 

PE  -  0.4,99 599 E  01  PSI 

SHOE  -C.MS3B8E-04  SLUGS/FT**3 
UE  ■  Q.4.62182E  04  FT/SEC 

li  -■  Imiimi  u,*y5/FT”J 

ME  -  0.1M745E  01 

NON-DIMENSIONAL  PRESSURE  COEFFTS.  NEWTONIAN-  1.173,  TANGENT  CONE-  1.249 
SHOCK  ANGLE  TO  SURFACE-  7.43  DEGREES 
REAL  GAS  VALUE,  CP  -0.1884E  01 

CALCULATION  OF  HEA1  TRANSFER  IN  ZERO  PRESSURE  GRADIENT 


REYNOLDS  ANALOGY  FACTOR  •  1.100 

REF.  .TEMP.,,  ECKERT  -  1464.1,  SOM-SHORT  -  1552.7,  DEG  K 
X  Q  BTU/SEC/FT**2  CFKS 


Q  BTU/$EC/F7**2 


TURBULENT 

TURBULENT 

LAMINAR 

O.SjOO 

Q.131718E 

03 

0.8.7637CE-02 

0.121495E 

03 

1.000 

0.111438E 

03 

0.7A1775E-02 

0.859104E 

02 

1.500 

0o 101586E 

03 

0.675890E-02 

0.7.01455E 

02 

2,000 

0.952803E 

02 

0.O33935E-02 

0.6.07478E 

92 

2.500 

0.907552E 

02 

0.6.09828E-02 

0.S.43345E 

02 

3.000 

0.&72758E 

02 

0.530677E-02 

0.496004E 

02 

3.500 

0.844764E 

02 

0.532052E-02 

0.4.59210E 

02 

4.000 

0.821512E 

02 

0.5A6582E-02 

0.429552E 

02 

CALCULATION  WITH  CORRECTION  FOR  CONI  Cl  TV 


X  IN 

RE 

PE 

CP(NEWT) 

CP(TC) 

Q(LAM) 

Q(S~S) 

Q(S-C) 

THETS 

O.SJOO 

0.1A38E 

05 

4.25 

1.152 

1.223 

120.081 

130.145 

73.086 

7.14 

1.000 

0.2O8SE 

05 

4.7.2 

1.131 
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83.8.90 

108.7,79 

64.908 

7.25 

1,500 

0.4344E 

05 

4.59 

1.110 

1.181 

67. 60S 

97,891 

60.699 

7,16 

2.Q00 

0,5^1'.E 

05 

4.40 

1.089 

1.159 

57,802 

90.6.1G 

57.7,24 

7.  (17 

2  ,  S".0 

C.72I/E 

05 

4.3,4 

1.069 

1.157 

51.131 

85.187 

55. 538 

6,9.8 

3.000 

0. 8.772  E 

05 

4.22 

1.050 

1.116 

46*093 

80.847 

53.820 

6.89 

5.300 

0.1026E 

06 

4.10 

1.030 

1.095 

42.136 

77.218 

52.165 

6,81 

4.000 

0.1177E 

06 

3.90 

1.011 

1.074 

38.9.13 

74.088 

51.109 

6.7,2 
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APPEIDIX  D 


Heat  Transfer  on  a  Sphere  and  a  Blunt  50°  Half-Angle  Cone 


Stagnation  point  heat  transfer 

The  theory  of  Pay  and  Riddell  (Ref.  1)  is  used  to 
predict  the  stagnation  point  heat  tranafer.  The  formula  used  is 


pp  o.i  du  2 f  1 

!  -  O.T63  Pr-»*«(p8y#)0.S(-a-ft)  (^-Sj #  [h#  -  hj 


Using  a  Modified  Sevtonian  theory,  the  Telocity  gra¬ 
dient  at  the  stagnation  point  of  a  body  vlth  radius,  Rq,  is 
given  by  : 

du  -  t2 

(«*> 

<!•  0  R*  /  9. 


where  p^  is  the  stagnation  pressure.  Subscript  e  refers  to 
conditions  just  outside  the  boundary  layer,  subscript  v  to  wall 
conditions  and  subscript  0  to  the  stagnation  point. 


The  total  enthalpy  h  was  assumed  e^ual  to  the  value 

C 

ia  the  nostle  reservoir,  calculated  for  real  nitrogen  free  the 
measurements  of  p0  and  Tq.  The  vail  enthalpy  hy  was  taken  with 
sufficient  accuracy  to  correspond  to  Tv  ■  2?3°S.  The  viscosity 
of  nitrogen  was  found  from  Sutherland's  formula  (Eon,  C.6)  with 
the  value  of  Tft,  the  static  temperature  behind  the  shock  calcu¬ 
lated  from  fc  taking  into  account  the  departure  from  a  perfect 
gas  due  to  vibrational  excitation  (Ref.  1*»).  The  stagnation 
point  pressure,  p.  ,  was  measured  and  the  densities  p  ,  u 
were  then  obtained  from  pt<  and  T#,  using  the  perfect  gi*e 
lav.  A  Prandtl  number  of  0.715  was  assumed. 


Curves  taken  from  the  Lees  similarity  theory  (Ref.  3) 
arc  used  to  predict  the  heat  transfer  rate,  in  relation  to  the 
stagnation  point  heat  transfer  downstream  of  the  stagnation 
point  for  these  body  shapes. 


. . . . . . . . . 


LONGSHOT  TEST  CONDITIONS,  RUN  NO.  285  AT  7,6.  IN.  .FROM  NOZZLE  EXIT, 
AND  4.6,  FROM  CENTRELINE 

T (MS)  PO(PSI)  TO(K)  PITOT(PSI) 

MACH  NO  POP(PSI)  TOP(K)  RE/FT 

PCPSI)  TOC)  RHO  V(FT/SEC) 

QD(LB/FT**2)  Q(BTU)  TT2R(K)  CONDENSATION 

0.0.00  0.585408E  05  0.240516E  04  0.806696E  01 

19.7,18  0.6.88585E  05  0.3.09855E  04  0.297866E  07 

0.158789E-01  0.3.93400E  02  0.181817E-94  0.827392E  04 

0.6.22343E  03  0.S87024E  02  0.264564E  04  0.3.59004E  02 


LONGSHOT  TEST  CONDITIONS,  RUN  NO.  291  AT  7.2.  IN.  .FROM  NOZZLE  EXIT 
AND  4.4.  FROM  CENTRELINE 

T  (MS )  PO(PSi)  TO(K)  PITOT  (PSD 

MACH  NO  POP (PS I )  TOP(K)  RE/FT 

PC  PS  I )  TOO  RHO  V(FT/S!  o 

QOCLB/FT**2)  Q(BTU)  TT2RCK)  CONDENS  JON 

0.0.09  0.3.76553E  05  0.222620E  04  0.5.36936E  01 

19.208  0.4.922 78E  05  Q.274521E  04  0.229039E  07 

0, llllSSE-Ol  0.3.67024E  02  0.136425E-04  0.7.78524E  04 

0.4.13437E  03  0.4.85692E  02  Q.231505E  04  0.346679E  02 


LONGSHOT  TEST  CONDITIONS,  RUN  NO.  292  AT  8.5,  III.  T  10M  NOZZLE  EXIT 
AND  4.9.  FROM  CENTRELINE 

T (MS)  PO(PSI)  TO(K)  PITOT(PSI) 

MACH  NO  POP(PSI)  TDP(K)  RE/FT 

P(PSI)  TOO  RHO  V(FT/SEC) 

00(LB/FT**2)  n(nTlj)  TT2R(K)  CONDENSATION 

0.0.90  0.579779E  05  0.2.34936E  04  0.292399E  02 

15.125  ft •  6.8 2 00 IE  05  0.392094E  04  0.8.13917E  C7 

0.9.76932E-01  *.&460C'._  02  0.681124E-04  0.813361E  04 

0.225301E  04  0.126923E  03  0.2563S5E  04  0.437353E  02 


TABLE  3  TYPICAL  TEST  .SECTION  CONDITIONS 
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TABLE  7  HEAT  TRANSFER  MEASUREMENTS  OH  “MODELS  D  & 
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